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INTRODUCTION 

(U)  This  report  describes  a computer  cataloging  technique 

that  is  capable  of  detecting  strong  signals.  It  was  developed 
for  use  in  analyzing  the  unprocessed  but  digitized  signal  data 
obtained  in  the  joint  US/UK  VERULAM /ANDREW  sea  trials.  Using 
the  technique  of  energy  averaging,  the  time  of  the  peak  arrival 
of  one  way  signals  and  bottom  reflections  are  detected  and 
recorded.  ^ 

(U)  The  methods  developed  in  this  report  supplement  previous 

cataloging  procedures^"  which  obtained  the  precursor  transmit  time 
directly  from  the  analog  tapes.  Rather  than  duplicating  previous 
work,  pertinent  data  is  inputted  into  the  computer  program.  From 
the  information  developed  by  these  methods  a master  catalog  of 
sonar  ping  cycles  can  be  obtained. 


TRACOR  Document  No.  SD/68-004-U,  "Summary  Report  - A Critique 
of  Analysis  Routines  for  VERULAM  Test,"  14  March  1968. 
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GENERAL  DESCRIPTION  OF  TECHNIQUE 


(C)  The  data  obtained  during  the  joint  US /UK  VERULAM  sonar 

tests  were  recorded  on  reels  of  12  channel  analog  tape.  In 
reducing  this  data  individual  channels  are  digitized  so  that 
digital  techniques  can  be  employed  to  further  process  the  data. 
One  of  the  tasks  of  highest  priority  is  the  processing  of  the 
two  way  received  signals  (echoes).  Because  of  the  inherent 
spreading  and  reflection  losses  of  the  transmitted  pulse,  the 
two  way  signal  has  a negative  signal- to-noise  ratio  prior  to 
signal  processing.  To  minimize  the  computer  time  required  to 
process  the  signal  data,  the  echo  arrival  time  is  estimated  by 
determining  the  travel  time  of  other  signals.  To  accomplish 
this  a computer  program  called  PEAK  TIME  was  written  to  permit 
searching  the  VERULAM  and  ANDREW  tapes  for  the  time  of  arrival 
of  certain  high  signal- to-noise  ratio  signals.  Knowing  the  time 
of  arrival  of  those  signals  would  permit  prediction  of  the  time 
of  arrival  of  low  signal- to-noise  echoes.  This  would  then 
minimize  the  duration  of  signal  correlation  analysis  necessary 
to  find  low  signal-to-noise  ratio  signals. 

(U)  The  computer  program  was  configured  to  search  for  the 

time  of  arrival  of  the  following  signals: 


VERULAM  Tape 


1.  Main  Pulse  Fathometer  Echo  Return 

2.  Main  Pulse  Slant-Bottom  Reverberation 

3.  Either:  Precursor  Direct-Path  Outbound /ANDREW 

Transponded  Direct-Path  Return 

Or:  Precursor  Bottom-Bounce  Outbound /ANDREW 

Iransponded  Direct-Path  Return 
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ANDREW  Tape 

1.  Precursor  One-Way  Direct  Path  Signal  (Optional) 

2.  Precursor  One-Way  Bottom- Bounce  Signal 

3.  Main  Pulse  One-Way  Direct-Path  Signal  (Optional) 

4.  Main  Pulse  One-Way  Bottom- Bounce  Signal 

COMPUTER  PROCESSING  METHOD 

(C)  The  original  VERULAM  and  ANDREW  data  tapes  were  analog 

recordings  of  the  unrectified  signal  received  at  each  ship.  These 
tapes  have  been  sampled  at  a frequency  of  6,250  Hertz  (.00016  sec. 
interval)  to  produce  digital  tapes  of  the  same  unrectified  signals. 
The  PEAK  TIME  program  first  adds  the  absolute  magnitude  of  a 
consecutive  group  of  samples  to  produce  an  output  proportional  to 
the  average  amplitude  of  the  unrectified  signal  over  the  duration 
of  the  sample  group.  A running  average  is  obtained  by  repeatedly 
adding  a new  sample  from  the  leading  edge  of  the  group  simultane- 
ous with  subtracting  the  oldest  sample  from  the  trailing  edge  of 
the  group.  In  searching  between  specified  limits  in  this  manner 
the  maximum  average  value  is  found  along  with  the  time  that  that 
maximum  occurred. 

(U)  To  search  the  entire  received  signal  data  for  each  ping 

cycle  is  considered  to  be  unnecessarily  time  consuming.  Instead, 
a calculation  is  made  to  determine  the  approximate  time  each 
signal  of  interest  is  due  to  arrive  and  then  to  confine  the 
search  to  a window  centered  about  that  time.  The  equations 
which  are  used  to  determine  the  window  locations  for  the  various 
signals  are  found  in  Appendix  A.  A derivation  is  given  for  the 
less  obvious  ones.  These  equations  are  used  to  locate  initially 
the  search  windows  in  a series  of  ping  cycles.  For  subsequent 
pings  the  windows  are  centered  about  the  actual  peak  arrival 
times  for  the  previous  ping.  If  the  peak  is  found  to  lie  outside 
of  a specified  percentage  width  of  the  basic  window,  as  measured 
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from  the  center  of  the  window,  the  peak  is  rejected  and  the  time 
previously  found  is  used.  Figure  1 illustrates  the  averaging 
interval  and  the  search  window  width  and  labels  many  values  with 
the  name  used  in  the  program. 

(U)  After  determining  the  time  of  each  peak,  the  program 

calculates  the  range  which  corresponds  to  that  time.  In  the 
case  of  the  fathometer  return,  the  range  is  the  water  depth. 

The  slant-bottom  range  is  the  distance  from  the  VERULAM  to  the 
bottom  along  the  assumed  straight-line  bottom-bounce  path  to  the 
ANDREW.  In  all  other  cases  the  range  calculated  is  the  horizontal 
range  from  the  VERULAM  to  the  ANDREW.  Average  range- rate  over 
the  most  recent  two  pings  is  also  calculated  for  the  VERULAM  to 
ANDREW  cases. 

(U)  A print-out  of  the  main  program  and  the  PEAKSK  sub- 

routine is  given  in  Appendix  B.  A complete  list  of  the  required 
program  inputs  by  Input  Card  Number  and  format  is  shown  in 
Appendix  C.  A reproduction  of  the  output  print-out  for  a typical 
run  is  given  in  Appendix  D along  with  an  explanation  of  the 
meaning  of  the  various  names  used. 

RESULTS 

(U)  The  program  was  test  run  using  VERULAM  digital  tape 

number  FT- 155  and  ANDREW  digital  tape  number  FT- 162.  These 
tapes  cover  the  LFM  pulses  from  ping  numbers  25  to  53  for  Run 
Number  684. 

(C)  Figure  2 shows  the  results  obtained  from  the  analysis. 

The  time  found  for  the  peak  signal  for  each  ping  is  shown  for 
the  fathometer,  the  slant-bottom,  and  the  bottom-bounce  trans- 
ponder/direct-path returns.  The  90  percent  window  width  limits 
within  which  the  signal  was  required  to  lie  are  shown  in  each 
case.  It  may  be  seen  that  the  fathometer  return  lies  well  within 
these  limits  in  every  case.  (The  fathometer  window  was  1.5  sec 
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VERULAM  DIGITAL  TAPE  NO.  FT- 155 


Window 

Width 

Limits 


PRECURSOR  BOTTOM-BOUNCE/TRANSPi 


DER  DIRECT-PATH  RETURN 


Window 

Width 

Limits 


SLANT- BOTTOM  RETURN 


FATHOMETER  RETURN 


Window 

Width 

Limits 


TRACKING  |OF  FATHOMETER,  S LANT- BOTTOM , AND  ANDREW 
TRANSPONDER  RETURNS  RECEIVED  BY  VERULAM 
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permitting  a deviation  of  .675  sec).  It  is  interesting  to  note, 
however,  that  the  fathometer  return  shifted  by  as  much  as  TL 
between  successive  returns  in  some  cases. 

(C)  The  slant-bottom  return  is  also  seen  to  lie  well 

within  its  prescribed  limits  of  .9  sec.  However,  this  is  not 
so  for  the  transponded  return  from  the  ANDREW.  Here  the  peak  was 
found  in  most  cases  to  lie  up  against  one  or  the  other  edges  of 
the  window,  indicating  that  there  was,  in  fact,  no  peak  within 
the  window.  A subsequent  examination  of  the  ANDREW  analog  tape 
indicates  that  the  transponder  delay  time  varied  in  an  erratic 
manner  and  this  explains  the  inability  of  the  program  to  find  a 
signal  peak  in  the  region  it  was  presumed  to  exist. 

(C)  Figure  3 shows  the  results  obtained  from  an  analysis 

of  the  ANDREW  tape  over  the  same  pings  as  covered  by  the  VERULAM 
tape  FT-155.  The  precursor  bottom-bounce  signal  was  required  to 
lie  within  very  narrow  limits  (+  .135  sec)  since  it  had  a pulse 
duration  of  only  0.1  sec.  It  is  seen  that  the  program  has  tracked 
the  signal  well  even  with  this  narrow  band.  The  main-pulse 
bottom-bounce  signal,  being  of  1.0  sec  duration  was  given  broader 
limits  (+  .9  sec).  The  program  appears  to  track  this  signal  also. 

CONCLUSIONS 

(C)  The  computer  program  appears  capable  of  tracking  all 

VERULAM  and  ANDREW  signals  well  except  those  transponded  by  the 
ANDREW.  This  failure  appears  to  be  caused  by  an  erratic  time 
delay  in  the  ANDREW  transponder.  Only  by  inputting  the  actual 
delay  time  for  each  ping  could  this  condition  be  rectified.  This 
would  require  examining  the  ANDREW  tape  to  determine  the  actual 
delay  time  for  each  ping  and  inputting  this  value  (ADELAY)  for 
each  ping  search.  The  program  would  have  to  be  modified  to  the 
extent  of  shifting  ADELAY  as  an  input  from  card  number  2 to  card 
numbers  5 and  beyond  as  necessitated  by  the  number  of  pings  being 
searched. 


CONFIDENTIAL  ROUGH  DRAFT 


ROUGH  DRAFT 


ANDREW  DIGITAL  TAPE  NOr  FT- 162 


^Window 

Width 

Limits 


MAIN  PULSE  BOTTOM- BOUNCE  SIGNAL 


Window 

Width 

Limits 


PRECURSOR  BOTTOM- BOUNCE  SIGNAL 


TRACKING  OF  PRECURSOR  AND  MAIN 
SIGNALS  RECEIVED  BY  ANDREW 


PULSE  BOTTOM- BOUNCE 


rougW  draft 


51  PING  NO 


UNCLASSIFIED  rough  draft 


msssf,.  2923  CANON  STREET,  SAN  OIEGO.  CALIFORNIA  92106 


APPENDIX  A 

TIME  AT  WHICH  SEARCH  WINDOWS  SHOULD  BE 
CENTERED  WHEN  LOOKING  FOR  SIGNAL  PEAKS  (U) 
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APPENDIX  A 

TIME  AT  WHICH  SEARCH  WINDOWS  SHOULD  BE 
CENTERED  WHEN  LOOKING  FOR  SIGNAL  PEAKS  (U) 

DEFINITION  OF  SYMBOLS 

R = VERULAM  to  ANDREW  horizontal  range  in  feet 
(Fig.  A-l) 

S = Straight-line  slant-bottom  distance  in  feet  along 
bottom-bounce  path  from  VERULAM  to  ANDREW 
(Fig.  A-l) 

h = Ocean  depth  in  feet  (Fig.  A-l) 

= Time  from  reference  to  beginning  of  transmission 
(Fig.  A- 2) 

tc  = Time  from  reference  to  center  of  transmission 
(Fig.  A- 2) 

ts  = Time  from  reference  to  end  of  transmission 
(Fig.  A- 2) 

t^  = ANDREW  transponder  delay  time  (Fig.  A-2) 

t^.  = ANDREW  transponder  pulse  duration  (Fig.  A-2) 

D 

C = Speed  of  sound  in  seawater  in  ft/sec 
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Ocean  Surface 


h (Ocean  Depth) 


Figure  A-l  DISTANCE  DIMENSIONS 


Pulse 

Transmission 


Pulse  Andrew 

Recognition  Transponder 

at  Andrew  Transmission 


Figure  A-2  TIME  DIMENSIONS 
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Time 
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A- 2 (U)  PRECURSOR  ONE-WAY  BOTTOM- BOUNCE  TIME 


*i  + h2 


2S  _ [r2  + Ah2 


One-way  travel  time  = =*  * ^ - — 

This  is  the  time  for  the  minimum  distance  path.  Multipaths 
will  cause  the  received  signal  to  reach  a maximum  at  the  end 
of  the  transmitted  signal.  Hence: 


R2  + 4h2 


LPC 1WBB 


If  the  travel  time  is  known,  the  range  can  be 
determined  by  solving  the  above  equation  for  R: 


[°2  0 


2 ai2 


PC1WBB  Te 


A- 3 (U)  SLANT-BOTTOM  REVERBERATION  TIME 

Due  to  the  vertical  width  of  the  beam,  the  bottom 
will  be  ensonified  at  different  ranges.  Therefore,  the 
reverberation  return  from  a small  time  increment  of  the  pulse 
will  be  spread  out  in  time.  Since  the  pulse  itself  has  a finite 
time  duration,  the  return  at  any  one  time  is  the  sum  of  the 
returns  from  different  parts  of  the  pulse  from  different  ranges. 
This  constitutes  a slice  through  the  "reverberation  ridge"  in 
Figure  A-3  at  the  time  in  question.  Such  a slice  is  shown  in 
Figure  A-4.  The  total  return  at  the  time  of  the  slice  will  be 
the  integral  under  the  curve.  However,  the  beginning  and 
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termination  of  the  pulse  may  truncate  the  curve  so  that  the 
integration  should  be  performed  between  limits.  The  location 
of  the  limits  will  slide  from  left  to  right  at  the  listening 
time  advances.  Shown  are  the  integration  limits  for  the 
receive  times  t^,  t2»  and  t^  marked  on  Figure  A- 3. 

(U)  Obviously,  the  maximum  will  be  achieved  when  the 

limits  straddle  the  curve  as  they  do  with  t^.  Time  t^,  then, 
is  when  the  bottom  reverberation  will  reach  its  maximum.  This 
time  is  equal  to  the  round-trip  travel  time  for  the  center  of 
the  beam  plus  half  the  pulse  length.  Thus: 


R2  + 4h2 


A-4  (U)  bottom-bounce  one-way  travel  time 
Same  as  Section  A-2  above: 


R2  + 4h2 


lBB1W 


■['*( 
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A- 5 (U)  PRECURSOR  ONE-WAY  DIRECT  PATH  TIME 


TPC1WDP  = § + t 
L e 


R ■ 0 [TPC1WDP  ■ Te 
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Figure  A- 3 SLANT- BOTTOM  REVERBERATION 

CHARACTERISTICS 


Figure  A- A SLANT-BOTTOM  REVERBERATION 

VERSUS  TIME 
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A- 6 (U)  DIRECT-PATH  ONE-WAY  TRAVEL  TIME 
Same  as  Section  A- 5 above: 

t „ R 4. 

^PIW  C Te 

R * c [TDP1W  ' Te] 

A- 7 (U)  FATHOMETER  ECHO  PEAK  TIME 

Analysis  is  similar  to  that  of  Section  A-3.  The  return 
from  the  bottom  resulting  from  a small  time  increment  of  pulse 
will  look  as  shown  in  Figure  A-5,  assuming  a narrow  beam  is  aimed 
directly  at  the  bottom. 

The  slice  through  the  bottom  echo  return  "ridge"  will, 
in  general,  look  as  shown  in  Figure  A-6.  If  the  receive  time  is 
relatively  early,  integration  limits  might  be  as  shown  as  t^ 
or  t2*  The  maximum,  obviously,  will  be  at  time  t^.  This 
time  is  equal  to  the  minimum  round  trip  travel  time  (for  the 
vertical  ray)  plus  the  time  from  the  reference  to  the  end  of 
the  pulse: 


T 


e 


A- 8 (U)  TRANSPONDER  TIME 

The  Transponder  Time  is  the  time  measured  from  the 
reference  for  the  pulse  to  be  transmitted,  travel  by  direct 
path  from  the  VERULAM  to  the  ANDREW,  be  transponded  by  the 
ANDREW  with  a time  delay  tA,  and  the  transponded  pulse  be 
returned  to  the  VERULAM.  The  direct-path  signal  outbound  will 
actuate  the  transponder  since  it  will  arrive  before  the  bottom- 
bounce  signal.  However,  the  transponded  signal  can  return  by 
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Figure  A- 5 BOTTOM  ECHO  FROM  SHORT 
DURATION  PULSE 
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both  direct-path  and  bottom-bounce  and,  thus,  there  will  be 
two  transponder  times.  In  both  cases  the  transponder  is 
assumed  to  be  triggered  by  the  leading  edge  of  the  transmitted 
pulse . 

(a)  Direct  Path  Out/Direct  Path  Return 

2R 

TXP  = TT  + Tb  CA  CB 


*-§[■ 


TXP  Tb  " tA  " CB 


(b)  Direct  Path  Out/Bottom-Bounce  Return 


TXP  = £ + Tb  + CA  + 


R2  + 4h2 


g2  [TXP  - Tb  - CA  - Cb]2  - 4h2 
2C  TXP  ‘ Tb  ' CA  " tB 


(c)  Transponder  Bottom-Bounce  Out/Direct  Path  Return 


R2  + 4h2 


+ Te  + CA  + Z + CB 


This  differs  from  (b)  only  in  that  Te  here  replaces 
t there.  In  (b)  it  was  assumed  that  the  transponder  was 
triggered  by  the  beginning  of  the  incoming  pulse.  Here  (from 
Section  A-4)  it  is  assumed  that  the  transponder  is  triggered  by 
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the  end  of  the  transmitted  pulse.  If  the  transponder  is 
triggered  a time  T after  the  arrival  of  the  wavefront  the 
equation  becomes: 


t . Isi+islL  + 1 + + ca  + cb 


Solving  for  range: 


c2  [TXP  - Tx  - CA  - 1 2 - 4h 

2C  [TXP  Tx  " CA  " Cb] 


(d)  Transponder  Bottom-Bounce  Out/Bottom  Bounce  Return 


2[R2  + 4h2l% 


]* 

— + <r  + t.  + t 


C2  ( 


x A B 


2 ,,2 


TXP  ' Tx  ‘ CA  ‘ CB  ) ' 4h 
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APPENDIX  B 

COMPUTER  PROGRAM  PEAK  TIME 
AND  SUBROUTINE  PEAKS K (U) 
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502  FORMA  I <lMl).5X.22HPRECORSflN  PULSE  LENGTH.  F10. ?.  «h  Sc  C . 

C21X.1/HMAIN  PULSE  LbNoTH , F 10 . 2 , 4H  SEC / / 1 4 X , 1 4 nPuL S t IWTcKVAL# 
CF 10  . ? » 4hi  StC*  23X,  15hA,vr>REW  OFF  T ! Mb  . F 1 0 . 2 * 4h  SfcC// 
CllX.l/hSAMPLI''JO  INTER VAL»F10.o*4H  ScC.?9X.9hWl'v0HHll)#Fl0.2. 
04 h Sp0//iox,or|WiNUw^(2).F10.?,4H  SEC*29X,yblwlNU8W(O),Fl0.2* 
C4n  SFU//49X#31HpEKCENl  WINDOW  WIDTp  ACCbPT  AbLt  * F 1U  • 2/// ) 
ShANfifc  ( l > sSQRT  l ( t .5*RANGFL>  **2  ) <►  < DEPTHL*nt:PTHL  > > 

FWUPpr  = (RANGFL/SflUNSP)-FPLbNrH(l)+PSTART(l)-ANDWFF 

pwDPMl=(RANGEL/SHUNSP)»PLcNrH(2)*PSrART(2)-ANUUFF 

fiWddpl'  = ((SRAN0fc(l)*2.0)/S0UNSP)+PLENTH(l)+PSTAKT(lJ-AND8FF 

RWbdMl=((SRANOE(l)*2.0)/S8UiMSP)fPLENTFi(2)*PSTART(2)-ANDdFF 

FATrtTM= (? . *DFPTHL/SOUNSP) ♦PLENTHC2) fPSTART  ( 2) 

SLANTF1=(  ( SKAN0b<  1 ) *2 . 0 ) /S8UNSP  ) ■*■ . 5*PLENTH  ( 2 ) -t-PSTAR  I ( 2) 

XUPOPl = ( ? . *R ANGEL/S0UNSP) +PSTART( 1 )+ADELAYf ADUKAT 

XUPOR  F = ( RANGFL/S0UNSP  J *®WddpT>ANDdFF-FPSTART(  1 J ♦ADt:LAY*AOUr(AT 

R ANuRl  ( 1 ) = R A N G b L. 


110 


502 


5o3 


11 

504 


RANG20 ( 1 ) = RANUc|_ 
A K A N m 1 ( 1 ) =RAN0cL 

AnANm211)sRaN0El 
xhi r ik=o . 

R R A i'  F 1 = 0 . 

ANOPP 1=0. 

ANDPm 1=0. 


R R A I M 1 = 0 . 

R c A 0 »OO,PlNGiM0,NVFIL!=,NAFlLE»SDlBd2,  IN.NLAST 

PH  I NT  -33o.PINuNo.NVFILc.NAFlLb 

PRINT  '5  31 

NOP  i P = 1 

JS AMP= J 

NfaSAMP=PLENTH(2)/SAMPNT 
pALF  TM=  . «3*W  in  JHW  ( 1 ) 

PRINT  o 3 2 

FciPtfA  i ( 1 H 0 » 6 H F ATH0M) 

CALL  PcAKSK  (FAJHI M.FATIMc.NVFILP) 

FRAiVGt(l^  = .5«60UNSP*(FATI'vlc-PSTART  (2>-PLENTH(2)  ) 
i»ALrTM=.H*WIND0w(2) 

PRINT  o 3 3 

F0R«aT(1H0*6HSL  HOT ) 


CALL  PcAKSK  (SLANlM.SLTIMt.NVFlLF) 

PrtI=5/.2957b*(ASIN(FAIIME/SLTiME)) 

NWSAkP=AOURAT/SamPNT 

Go  IP  (li,12).SDldB2 

PRINT  p34 

forma i (ih0»6hx  dpop) 

CALL  PcAKSK  ( XDPDPT , XPl T I M . NvF I LE ) 
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R A N lie  1 ( £ ) sRANlitK  i ) 

RA.\Uipi(l)  = ( .5*b*UINSP)*(XPiriM-PSTART(l) -ADEL AY-ADO* AT) 
RKAIcl=(RANGFl(l)-RANGEl(2)>/<PULSNT*Fl.@AT(lN>> 

Go  1 P)  It 
12  PRINT  •>  0 5 

505  F t'RMAl  < lHO.orlX  BRDP) 

CAll  ^tAKSK  ( xuprioT,  XP2T  I M»  NVF  ILE ) 

RANG?0(2)sRaNG23<1) 

tip.  Ilo  = XP^TIM-PSTART<1  ) - AJCLAY-ADURAT 

RANu?o(l>  = < < (a0lJNSP*TlMTW0)**2)-{4.O*FRANGE<l>*FRANG6<l)))/ 

C< 2 .0*bbUNSP*T  IMTWio) 

RRA  T?o  - ( RANG20(  1 ) -RANG23I  2)  )/(  PUl.SNT*FLbAT(  IN>  ) 

14  KsriiP  = 2 
jaAMp=0 

NcisAMR  = PLcNTN(l)/SAMPNT 
HALFTM=.S*WINJ0W(O) 

Gb  T p ( lx , 17 ) » SD1O02 

16  PRINT  26* 

5 0 o FaiRMAl(lH0»6hbR  P l)  P ) 

CALL  PtAKSK  (bRDPPT#  A.vDPPT*NAFlLF) 
a N 0 R P 1 = A\iDPPT*A\nciFF 

17  PRINT  P07 

507  FbHNAl(lHO*6H0W  Pori) 

CALL  RcAkSK  ( ri  W R R P T , aNHBPT  » N A F I L F ) 

A N a d d I = AnRbPT*AND«IFF 

NfflbAMP  = PLt:NTrt(2)/bAMPNT 
f.ti  Ip  1 1 ? > i 9 > * b D 1 to  9 2 

18  PRINT  ^OP 

508  F0AMaH1HO»6Hciw  M D P ) 

CALL  PcAKSK  (0WDPMT,ANnPMT*NAFlLE) 

ANDRm I = A \ D P M T * A N D b F F 
ARANMl  ( 2 ) = AR  A.MMi  ( 1 ) 

ARANm1(1)=S0UNSP»iANDPMT-PPTART(?;-PLENTH(2)) 

Pr  A Im1=(aPANMlI1)-ARANM1(2) )/(Pulsnt*flsaT( in) ) 

19  hALF  T«=.5*RI,\D0W<  1) 

PrINT  o6q 

509  FbKMAl(lhU»6H0W  M b B ) 

call  peaksk  (wwaPMTf aN3bmt»naf Ilf) 

ANbDMf=ANBriMT‘FA|\iUbFF 

AKANm2(2)sARANM2(1) 

,C.  5^  ^ 
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1 


AHAN>'C(l)=S0F?T(((bWUNbP*(ANd«MT-PbTART<2)-PLP>*Tr1(^>))**2)- 
C(4.*(hKANGt(l)**/*J  ) ) 

RRA  Hl2s( aRANm21 1 > -AKANM2( 2>)/ ( PULSN  F*FL0A  I ( IN)  ) 

PK 1 NT  340 

PnlNT  370,PSTART<l>,P5TARTC2>»FATlMb,FPANGE<l)»SLTiME*Prtl. 

CXH1 1 I M. RaNGEI ( 1 ) , bRATcl 
PKlNT  341 

prclNT  3/l,XP2TlM,KANG23(l>»rtRAT23»ANnPPT,ANB0Pl*ANupMT. 
CAKANMl(l).HRATMl,ANbBMT»AKA(MM2(l)*H^ATM2 
5,50  FidHMAf  (1H1*//»12H  PING  NR.  = * I 3 . 10  X , 1 9HV  £ R AL  w M F I u&  N0.  = .115. 
ClGX.lOrlANDREW  F I L t:  NFI.  = >16//) 

b6 1 FtdNMA  H IhG  » 10X  . fcHNRSArtP  ,5X.  6H!>AMPl  M » 4X  , &HJSTAH  T , 6X#  6rtST  AK  I M»  5*  » 

C5HJS>ToP,(SX.6HST0PTM.9X,6HSlGTIM,lOX.4HJMAA,9X.oriPtAKTM,/) 

540  FGHMa  I (1hO»//X/  / , 7X.9HPSTAHI  (1)  »?X*9HPSTahT  (?)  ,4X,orlFATlMe,3X* 
CftnFr<A'NGfc,5X,AhSLTiMfc,6X.3r(Pm.6X,6HXPlTlM.5X.6hKANutl,5X.0hHRATEi; 

541  F0HMAl(lHO./.6HxP^TIM.5X.6FIHANG23.5X,6HRRAT23.4X,6f1ANiJPPT.4X. 
C6nANRdPT,4X,6HANDPMT.3X,6r1ArtANMl,5X,6HRRATMl,«fX,6rl*NddMT»5X. 
C6rtAriAi'Hvl2, 5X,  6HRRA  I M2) 

560  FpHMaT  (lhO»4F20.6) 

570  F0RMa I (lH0»Fl2.2*Fll.2,F12*2»Fl3.2»Fl0.2»F10.2#Fli*2*F12.2» 

CF1U . ? » / ) 

571  F0RMa1v1hG»F6.2,F12.?,F1O.2,F1O.2#F1O.2,F1O.2*F12.2,F1O.2» 
CF10.?*Fll.2.F11.2#////> 

FATnTM=FATlME 
SLANTS s5|.  T I ME 
xupup i =xpi r i m 
XOPdp I = X P 2 T I M 
PirtUPp  T = ANDPPT- AND0FF 
(.  ivrirtP  i = ANibdPT"A(\|DbFF 
p m L)  P M I = AN'DPMT"Ai<inbFF 
(•)  a a d M 1 = ANdBMT'ANDaFF 

IF ( NLAbT ) 110.110,120 

120  CBNT  J N U fc 
END 


at* 
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SUbrtdUriNP  pfcARbK  ( S I G I I M » PE  AK  T M , Nl ) 
n 1 M c N 0 I 0 M I S I 0 ( 6 3 0 0 ) 

C o M M p I\  hALFTM.HEWJNT.SAMPNI  »Nc)SAMP»NSFiIP,JSAMP 
IF  ( b I G T I «• ) 30,81.81 

M JSTAP  1 = ( ( SliiT  IM-hALFTM  ) / SAMPNl  ) - < .5*<  FLWAT  ( NOSAMP)  > ) 
I K ( JS I APT ) 82,82,03 

02  JOT  AH  I =1 

03  jSTbP=(  (sICiTIM-fh4LFTH)/SAMPNTJ-f(  . 5# ( FL0AT ( N0S AMP ) ) > 

S'  I A H T ivi  = b A H P N T * ^ F L ft  A T ( J S T A H I ) ) 
SlbPTM=SAMPNT*tFLbAT(JSTSPi) 

SAMPTivl  = bAMPNT*(FLB'AT(NflSAF'P>) 

i r w t a l = o 

IF ( JS i SP- JSAMP)  O0.80.fa6 


JSTAH  f = Jfa  AMP 


06  I F ( JS I AKT  . LT . JfaAMP  ) JSTArtF; 
jcNLir  JST  ART-JfaAMP 
Go  Iftlo,6,50,o0>,i\SF1lP 

5 CALL  PdST(N) 
r a l l olkpt 

\ S FI  ip  = 0 

50  * = u 

31  Call  KucTPT  ( Si  1>V  AL  , NDK  I LE  ) 

J fa  AMP  = JS  AM  P + 1 
K = K + 1 

I F ( IM H F I LF  ) 

60  IF  (A-JcNFi) 

6 CALL  p 0 S T M 3 ( N ) 

PALL  ULPk3 

\ fa  ,i  i P = 4 
oO  K = 0 

36  CALL  KobTM3(SlGVAL.NDFlLF J 
jsAMp= JSamP+1 
K = K + j 

]F  (MFiF  ILF) 

6 5 IF (6-JtNn) 

9 n » 2 iM  S = i , N 0 S A M p 

GO  IF  (41,42. 41, 4 2 ) » N S H I P 
*tl  C-AuL  KiibTPKSlUVAL.NDFlLE) 
pa  i p 43 

42  (ALL  KotTM3(SI0VAL,NDFlLE) 

43  JfaAMP= JS AMP+1 
I F l MI)F  i LF  ) 

10  lSlb(NS)rAHS(faiGVAL) 

I 1a1ai-=ITKTAL*ISIu(NS) 

2 c b N 1 I iM  U E 

IPK  I ft  I S I Tfa'f  At. 
jMAXrJSAMP" ( NbfaAMP/2  > 
PcAKTH=(FL0AT<JMAX>)*SaMPNT 
15  D b 4 iMS  = i,NdSAMp 


30.30,80 

31,9.9 


35,35,60 

36,9,9 


10,10,80 
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I n-)UL  = imAl-lStLi(NS) 
pci  I R (45,46,4t>,4o),NS>HIP 

45  CALl  KiitTPT  (SlbVAL.NDK  ILE) 

Gd  Tfl  47 

46  CALL  KuETM3(SIGv  AL.NDF  ILF. ) 

47  jSAMps JSAMP*1 

I f ( Nllt-  ILF) 

20  IbIli(NS)aAe<S(i>IoVAL) 

I I el  i al=  i twTal*  l $ ! i><  NS ) 

I F ( ITiafAL-lPK  reij) 

25  IPM  p I = UdTAL 
JMAXsu6AMP-(NUSaMP/2) 

26  IP  ( JSAMP-JSTBK)  4»  75»  7t> 

4 C W N T I N U fc 

GO  Tp  15 

75  Pca6TM= i FLOAT ( JmA A ) ) *SAMPNT 

PH  I NT  106  . NO  S' amp,  6AMPT.H,  jS  lART,STARTM,JSTfe!P,STciPTrt»SIi»TlM,JMAX, 
0 P t A K T rt 

IF  ( ( AOSl S I GT I M-PEAKTM  > > - < PtHCNT*HALFTM  > ) 60*  PO, 

76  PtAMM  = SlGTlM 

PH  1 NT  104  » PE  A KTm 
CO  CONTINUE 

04  FoHMaT(1h0#10C»X,F1?.4,/) 

05  FOHMaT  (1H-F.9H  SlGNAL..,I6#H2.4,I10»F12.4,U0»F'12.4.Fl5.4,I15» 
CF14  . a ) 

Re  TURN 
END 
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APPENDIX  C 


PROGRAM  INPUTS  (U) 

C-I  (U)  INPUT  CARD  NO.  1 (FORMAT  8110) 

NORUN  = VERULAM/ANDREW  Run  Number.  Columns  1-10. 

NVTAPE  = VERULAM  Digital  Tape  Number.  Columns  11-20. 

NATAPE  = ANDREW  Digital  Tape  Number.  Columns  21-30. 

NOFILE  = Number  of  Files  Being  Searched  in  This 
Computer  Run.  Columns  31-40. 

C-2  (U)  INPUT  CARD  NO.  2 (FORMAT  8F10.4) 

VS PEED  = VERULAM  Speed  in  Knots.  Columns  1-10. 

VDEPTH  = Depth  of  VERULAM  Transducer  in  Feet.  Columns 
11-20. 

VDEPRS  = VERULAM  Transducer  Center-Beam  Depression 
Angle  in  Degrees.  Columns  21-30 

DEPTHL  = Logged  Value  of  Ocean  Depth  in  Feet.  Columns 
31-40.- 

RANGEL  = Logged  Value  of  VERULAM- AN DREW  Range  in  Feet. 
Columns  41-50. 

APSECT  = ANDREW  Aspect  Angle  in  Degrees  as  Viewed 
From  VERULAM.  Columns  51-60. 

ADEPTH  = ANDREW  Depth  in  Feet.  Columns  61-70. 

ADELAY  = ANDREW  Transponder  Delay  Time  in  Seconds. 
Columns  71-80. 

C-3  (U)  INPUT  CARD  NO.  3 (FORMAT  8F10.4) 

ADURAT  ■ Duration  in  Seconds  of  ANDREW  Transponded 
Pulse.  Columns  1-10. 
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PLENTH(L)  - VERULAM  Precursor  Pulse  Length  in  Seconds. 
Columns  11-20. 


PLENTH(2)  - VERULAM  Main  Pulse  Length  in  Seconds. 
Columns  21-30. 

PSTART(l)  = VERULAM  Precursor  Start  Time  in  Seconds 
Measured  From  the  Time  of  the  First 
Sample  in  the  File.  Columns  31-40. 

PSTART(2)  = VERULAM  Main  Pulse  Start  Time  in  Seconds 

Measured  From  the  Time  of  the  First  Sample 
in  the  File.  Columns  41-50. 


SAMPNT  = Interval  in  Seconds  Between  Successive 
Digital  Samples.  Columns  51-60. 

PULSNT  = Interval  Between  Successive  VERULAM  Pulse 
Transmissions  in  Seconds.  Columns  61-70. 


ANDOFF  = The  Length  of  Time  in  Seconds  Between  the 
First  Sample  in  VERULAM  File  and  the  First 
Sample  in  the  ANDREW  File.  Columns  71-80. 


C-4  (U)  INPUT  CARD  NO.  4 (FORMAT  8F10.4) 

XLEVEL  = VERULAM  Transmit  Center-Beam  Source  Level. 
Columns  1-10. 


1 


BWIDTH  = Receive  Bandwidth  in  Hertz.  Columns  11-20. 

PERCNT  * Percent  of  Search  Window  Width  Located  in 
Center  of  Window  Within  Which  Peak  Value 
Must  Lie  to  be  Accepted  as  Peak  Time. 

SOUNSP  = Speed  of  Sound  in  Seawater  Measured  in 
Feet  per  Second.  Columns  31-50. 

WINDOW (1)  = Search  Window  Width  in  Seconds  Used  in 
Searching  for  VERULAM  Fathometer  Return 
and  ANDREW  Bottom-Bounce  Main  Pulse  Signal. 
Columns  41-50. 
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WINDOW (2)  « Search  Window  Width  in  Seconds  Used  in 
Searching  the  VERULAM  Tape  for  the 
Following  Returns: 

1.  Slant  Bottom  Reverberation 

2.  Direct  Path  Outbound/ANDREW  Transponded 
Direct-Path  Return 

3.  Bottom-Bounce  Outbound  ANDREW  Transponded 
Direct-Path  Return 

Columns  51-60 

WINDOW (3)  ■ Search  Window  Width  in  Seconds  Used  in 

Searching  the  ANDREW  Tape  for  the  Following 
Signals : 

1.  Precursor  Direct-Path  Signal 

2.  Precursor  Bottom-Bounce  Signal 

3.  Main  Pulse  Direct-Path  Signal 
Columns  61-70. 

C-5  (U)  INPUT  CARD  NO.  5 AND  BEYOND  (FORMAT  8110) 

(A  card  for  every  VERULAM /ANDREW  file  to  be  searched) 

PINGNO  « VERULAM  Ping  Number  Being  Investigated. 
Columns  1-10. 

NVFILE  = File  Number  on  VERULAM  Tape  Being  Used 

Which  Corresponds  to  PINGNO.  Columns  11-20. 

NAFILE  = File  Number  on  ANDREW  Tape  Being  Used  Which 
Corresponds  to  PINGNO.  Columns  21-30. 

SD1BB2  ■ 1 for  Surface-Duct  Operation.  (Beam 

depression  angle  at  or  near  zero  degrees.) 

■ 2 for  Bottom-Bounce  Operation.  (Beam 
depressed  well  away  from  horizontal.) 

Column  40. 
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IN 


NLAST 


- Number  of  Transmitted  Pulses  Since  Last 
File  Searched.  Columns  41-50. 


■ On  the  Input  Card  for  the  last  VERULAM/ 
ANDREW  file  to  be  searched  an  integer  i 
should  be  placed  in  Column  60  to  cause 
the  run  to  terminate. 
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APPENDIX  D 
PROGRAM  OUTPUTS  (U) 

(U)  The  program  initially  prints  out  the  following  basic 

information  on  the  run: 

RUN  NUMBER  - The  VERU LAM /ANDREW  Run  Number 

NUMBER  OF  FILES  - The  number  of  pings  being  analyzed 

in  this  computer  run 

VERULAM  DIGITAL  TAPE  NUMBER 

ANDREW  DIGITAL  TAPE  NUMBER 

VERULAM  SPEED  IN  KNOTS 

ANDREW  ASPECT  with  respect  to  the  VERULAM  in  degrees 

VERULAM  SOURCE  DEPTH  in  feet 

ANDREW  DEPTH  in  feet 

VERULAM  DEPRESSION  ANGLE  in  degrees 

TRANSPONDER  DELAY  TIME  in  seconds 

INITIAL  OCEAN  DEPTH  in  feet 

TRANSPONDER  PULSE  DURATION  in  seconds 

TARGET  RANGE  - VERULAM  to  ANDREW  range  in  feet 

TRANSPONDER  LEVEL  in  dB 

SOUND  SPEED  - Assumed  speed  of  sound  in  seawater 
in  feet/sec 

BANDWIDTH  in  Hertz 

PRECURSOR  PULSE  LENGTH  in  seconds 

MAIN  PULSE  LENGTH  in  seconds 

PULSE  INTERVAL  in  seconds 

ANDREW  OFFTIME  - Time  in  seconds  after  VERULAM  reference 

time  before  ANDREW  recorder  was  turned  on 
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SAMPLING  INTERVAL  - Time  in  seconds  between  consecutive  samples. 

WINDOW (1)  - Time  interval  in  seconds  over  which  the  following 
signal  searches  will  be  made: 

VERULAM  Tape:  Main  Pulse  Fathometer  Echo 

ANDREW  Tape  : Main  Pulse  Bottom-Bounce  Signal 

WINDOW (2)  - Time  interval  in  seconds  over  which  the  following 
signal  searches  will  be  made: 

VERULAM  Tape:  Main  Pulse  Slant-Bottom  Reverberation 

Precursor  Direct-Path  Transponder  - 
Actuated  Direct  Path  Return 

Precursor  Bottom-Bounce  Transponder  - 
Actuated  Direct  Path  Return 

WINDOW (3)  - Time  interval  in  seconds  over  which  the  following 
signal  searches  will  be  made: 

ANDREW  Tape:  Precursor  Direct-Path  Signal 

Precursor  Bottom-Bounce  Signal 

Main  pulse  Direct-Path  Signal 

PERCNT  - The  percent  of  the  full  window  width  within  which  the 

peak  average  signal  value  must  lie  to  be  acceptable.  See 
Figure  1.  If  peak  lies  within  this  limit,  it  is  used  as 
center  of  window  for  corresponding  search  on  next  ping. 

If  not,  previously  established  window  center  is  used. 

(U)  Following  the  listing  of  the  basic  information  identifying 

the  run  the  program  prints  out  solution  values  found  for  each  ping 
transmission  within  the  run.  These  are  identified  by  the  PING  NO. 
as  well  as  VERULAM  FILE  NO.  and  ANDREW  FILE  NO.  associated  with 
the  ping.  Immediately  beneath  these  headings  nine  calculated 
values  are  printed  under  the  following  names: 
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NOS  AMP 


SAMPTM 


JSTART 


STARTM 


JSTOP 


The  number  of  samples  (at  .00016  sec  intervals) 
over  which  a running  absolute  amplitude  average 
of  the  unrectified  signal  is  maintained. 

The  time  duration  in  seconds  represented  by 
NOS AMP. 

The  sample  number  for  the  first  sample  in  the 
first  NOSAMP  averaging  interval  for  a window 
search. 

The  time  in  seconds  at  JSTART.  Time  is  measured 
from  the  reference  which  is  established  by 
PSTART(l) , the  time  of  commencement  of  the 
precursor  transmission. 

The  sample  number  for  the  last  sample  in  the 
last  NOSAMP  averaging  interval  for  a window 
search . 


STOPTM  = The  time  in  seconds  at  JSTOP. 


SIGTIM 


JMAX 


PEAKTM 


The  time  about  which  the  window  search  was 
carried  out. 

The  sample  number  for  the  middle  sample  of  the 
NOSAMP  interval  which  had  the  highest  average 
amplitude  value  over  the  particular  window 
search. 

The  time  at  JMAX.  If  this  time  deviates  from 
SIGTIM  by  more  than  the  inputted  percentage 
value  (PERCNT)  of  the  half-window  width,  the 
calculated  value  of  PEAKTM  is  rejected.  It  is 
replaced  by  SIGTIM.  When  this  occurs  it  is 
shown  by  the  repeat  print-out  of  the  SIGTIM  below 
and  slightly  to  the  right  of  the  PEAKTM  print- 


UNCLASSIFIED  rough  draft 


UNCLASSIFIED  rough  draft 


2923  CANON  STREET.  SAN  DIEGO,  CALIFORNIA  92106 


(U)  Window  searches  may  be  conducted  for  any  of  eight 

different  signal  returns.  These  are  identified  at  the  left  under 
the  following  names: 

FATHOM  - The  fathometer  or  vertical  echo  return 
from  the  bottom  to  the  VERULAM  from  the 
main  pulse 

SL  BOT  - The  reverberation  return  from  the  bottom 

to  the  VERULAM  from  the  center-beam  of  the 
main  pulse 

X DPDP  - The  precursor  direct-path  outbound- transponder 
direct-path  return-signal  to  the  VERULAM 

X BBDP  - The  precursor  bottom-bounce  outbound- 

transponder  direct-path  returned  signal 
to  the  VERULAM 

OW  PDP  - The  one-way  VERULAM  to  ANDREW  precursor 

direct-path  signal 

OW  PBB  - The  one-way  VERULAM  to  ANDREW  precursor 
bottom-bounce  signal 

OW  MDP  - The  one-way  VERULAM  to  ANDREW  main-pulse 
direct-path  signal 

OW  MBB  - The  one-way  VERULAM  to  ANDREW  main  pulse 

bottom-bounce  signal 

(U)  The  results  of  the  search  are  set  forth  in  the  last 

part  of  the  print-out  for  each  file.  The  names  here  indicate: 

PSTART(l)  - The  time  of  commencement  of  the  precursor 
pulse,  measured  in  seconds  from  the  first 
sample  on  the  VERULAM  file.  Since  the 
start  of  the  precursor  is  usually  the 
reference  time  from  which  all  other  times 
are  measured,  PSTART(l)  is  usually  zero. 
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PSTART(2)  - 


The  time  in  seconds  from  the  first  sample 
on  the  VERULAM  file  at  which  the  main  pulse 


commences 


FATIME  - The  time  in  seconds  at  which  the  fathometer 
return  signal  peak  was  found  to  occur 

FRANGE  - The  ocean  depth  in  feet  corresponding  to 
FATIME  at  the  sound  speed  used 

SLTIME  - The  time  in  seconds  at  which  the  slant- 

bottom  reverberation  peak  was  found  to  occur 

PHI  - The  straight  line  depression  angle  resulting 

from  FATIME  and  SLTIME 

XPlTIM  - The  time  in  seconds  at  which  the  precursor 
direct-path/direct  path  transponder  signal 
peak  was  found  to  occur 

RANGE1  - The  horizontal  range  corresponding  to  XPlTIM 

RRATE1  - The  range-rate  in  ft/sec  resulting  from  the 

last  two  calculated  values  of  RANGEl. 

(First  calculation  uses  RANGEL  and  RANGEl.) 

XP2TIM  - The  time  in  seconds  at  which  the  precursor 
direc t-path/bottom-bounce  transponder 
signal  was  found  to  occur 

RANG23  - The  horizontal  range  in  feet  corresponding 
to  XP2TIM 

RRAT23  - The  range-rate  in  ft/sec  resulting  from  the 
last  two  calculated  values  of  RANG23. 

(First  calculation  uses  RANGEL  and  RANG23.) 

ANDPPT  - Time  of  arrival  in  seconds  at  which  the 

direct-path  precursor  signal  peak  was  found 
to  arrive  at  the  ANDREW 
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ANBBPT  - Time  of  arrival  in  seconds  at  which  the  bottom- 
bounce  precursor  signal  peak  was  found  to 
arrive  at  the  ANDREW 

ANDPMT  - Time  of  arrival  in  seconds  at  which  the 
direct-path  main  pulse  peak  was  found  to 
arrive  at  the  ANDREW 

ARANMl  - The  horizontal  range  in  feet  corresponding  to 
ANDPMT 

RRATMl  - The  range-rate  in  ft/sec  resulting  from  the 

last  two  calculated  values  of  ARANMl.  (First 
calculation  uses  RANGEL  and  ARANMl.) 

ANBBMT  - Time  of  arrival  in  seconds  at  which  the 

bottom-bounce  main  pulse  was  found  to  arrive 
at  the  ANDREW 

ARANM2  - The  horizontal  range  in  feet  corresponding 
to  ANBBMT 

RRATM2  - The  range-rate  in  ft/sec  resulting  from  the 

last  two  calculated  values  of  ARANM2 . (First 
calculation  uses  RANGEL  and  ARANM2 . ) 
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